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基于 核 函数 与 马 氏 距离 的 FCM 图 像 分 割 算法 


王 ” 燕 ， 才 祥 惠 ， 段 亚 西 
(兰州 理工 大 学 计算 机 与 通信 学 院 ， 兰 州 730050) 


摘 要 : 针对 模糊 聚 类 算法 邻 域 信息 与 空间 信息 利用 率 低 易 受 噪声 影响 的 问题 ， 提 出 一 种 结合 核 隙 数 与 马 氏 距离 的 
FCM 算法 ， 即 FCMKM 算法 。 首 先 ， 将 图 像 像 素 点 由 低 维 空间 通过 核 范 数 非 线性 映射 到 高 维 空间; 然后 ， 利 用 马 
氏 距 离 替换 原 有 的 欧式 距离 作为 高 维 空间 距离 量度 ; 最 后 ， 利 用 改进 后 的 算法 对 图 像 进行 分 割 。 为 验证 FCMKM 算 
法 的 性 能 ， 选 取 Bezdek 划分 系数 、Xie_Beni 系数 、 重 构 错 误 率 、 运 行 时 间 、 迭 代 次 数 五 个 评测 指标 作为 对 比 实验 
的 评价 标准 。 实 验 结果 表明 , 与 传统 FCM 算法 、 基 于 核子 数 的 FCM 算法 、 基 于 马 氏 距离 的 FCM 算法 相 比 , FCMKM 
算法 能 有 效 的 提高 模糊 聚 类 算法 的 抗 噪 性 。 
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Image segmentation of FCM algorithm based on kernel function and Markov distance 


Wang Yan, Qi Xianhui, Duan Yaxi 
(College of Computer & Communication, LanZhou University of Technology, Lanzhou 730050, China) 


Abstract: In order to solve the problem that the low utilization rate of the neighborhood information and spatial 
information led to vulnerability to noise of fuzzy clustering algorithm, This paper proposed a fuzzy C-means algorithm 
combining kernel function and Mahalanobis distance(FCMKM) . Firstly, it non-linearly mapped the image pixels from 
low-dimensional space to high-dimensional space through the kernel function. Then, it replaced the original Euclidean 
distance with Mahalanobis distance as a high-dimensional spatial distance measurement. Finally, it used the improved 
algorithm to segment the image. The paper selected five evaluation indexes Bezdek partition coefficient, Xie_Beni 
coefficient, reconstruction error rate, running time and iteration number as evaluation criteria of comparative experiments to 
verify the performance of the FCMKM algorithm. Experimental results show that compared with traditional FCM algorithm, 
kernel function based FCM algorithm and markov distance based FCM algorithm, FCMKM algorithm can effectively 
improve the anti-noise performance of fuzzy clustering algorithm. 
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0 ”引言 数 的 特性 ， 不 能 很 好 地 解决 边界 噪声 问题 ， 使 得 到 的 分 割 结 
六 果 存 在 误差 。 文 献 [7] 通 过 将 像素 相关 性 引入 到 模糊 因子 中 ， 
到 图 像 分 割 在 图 像 分 析 的 预 处 理 阶段 中 具有 十 分 重要 的 作 ”提高 了 模糊 聚 类 算法 抗 噪 性 ， 但 是 处 理 效率 较 低 、 算 法 时 间 
用 ， 是 图 像 后 期 分 析 的 基础 。 图 像 分 割 是 根据 图 像 的 各 方 男 复杂 度 较 高 。 文 献 [8] 引 入 马尔 可 夫 随 机 场 函 数 ， 能 有 效 地 提 
属性 将 图 像 分 为 多 个 互 不 目 具 有 相同 属性 的 区 域 ， 各 区 ”高 对 模拟 噪声 的 抗 噪 能 力 ， 但 是 对 处 理 含有 自然 噪声 的 图 像 
域 压 加 为 整个 图 像 ， 区 域 盟 性 相似 性 高 ， 区 域 之 间 差 异 ”时 提高 不 明显 。 文 献 [9 通过 引入 邻 域 相 对 粒 来 计算 隶属 度 ， 
性 明显 外。 该 算法 提高 了 模糊 聚 类 算法 的 抗 噪 性 , 但 是 增加 了 运算 时 间 ， 
模糊 聚 类 算法 是 模糊 理论 与 聚 类 算法 结合 产生 的 新 算 ”使 算法 的 整体 性 能 下 降 。 文 献 [10] 采 用 基于 广义 高 斯 分 布 的 
法 , 通过 对 数据 集 进 行 模糊 划分 , 得 到 合适 的 聚 类 结果 .Dunn 模糊 聚 类 算法 对 模糊 聚 类 结果 进行 边界 修整 ， 但 
等 提出 模糊 C 均值 算法 (FCM)， 并 且 Bezdek 等 人 BI 首次 将 ”中 的 边界 与 实际 边界 存在 误差 ， 分 割 过 程 丢 失 边 界 信息 ， 仍 
模糊 聚 类 算法 应 用 到 图 像 分 割 领域 。 近 年 来 ， 许 多 学 者 在 基 ”无 法 解决 噪声 问题 。Hussain 等 人 00 提 出 一 种 基 
于 模糊 聚 类 的 图 像 分 割 方面 作出 了 大 量 的 成 果 。 为 提高 模 煌 部 信息 的 FCM (KWFLICM) 算法 ， 该 算法 利用 核 函 数 能 有 
聚 类 算法 的 抗 噪 性 ， 文 献 [外 通过 引入 邻 域 隶属 函数 ， 加 强 了 ” 效 地 提高 模糊 聚 类 算法 的 抗 噪 性 ， 但 是 在 处 理 自然 图 像 时 存 
邻 域 对 像素 隶属 度 确 定 的 影响 ， 但 应 用 于 自然 图 像 和 医学 图 。 ”在 分 割 不 完整 现象 。 Zhao 等 人 (3 提出 一 种 结合 马 氏 距离 的 
像 分 割 时 却 出 现 了 丢失 空间 信息 的 问题 。 文 献 [3] 通 过 引入 软 。 FCM (FCM-M) 算法 , 该 算法 能 有 效 地 提高 离散 数据 的 聚 类 
集 概 念 ， 得 到 改进 的 软 集 模糊 聚 类 分 割 算法 ， 但 由 于 CT 图 “效果 ， 提 高 了 算法 的 抗 噪 性 ， 但 是 在 应 用 于 图 像 分 割 时 ， 不 
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像 特点 ， 使 得 实验 结果 图 像 边 界 噪声 明显 、 抗 噪 性 能 较 差 。 能 有 效 地 分 割 目标 区 域 。 
为 解决 噪声 对 边界 分 割 的 影响 ,文献 [6] 提 出 了 一 种 多 目标 进 为 了 解决 模糊 聚 类 算法 噪声 问题 ， 本 文 提 出 一 种 结合 核 
化 模糊 聚 类 算法 ， 由 于 模糊 聚 类 算法 具有 过 度 依 赖 隶 属 度 函 函数 与 马 氏 距离 的 模糊 聚 类 算法 (fuzzy C-means algorithm 
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combining the kernel function with the Mahalanobis distance, 
FCMKM)， 即 利用 核 函数 与 马 氏 距离 结合 改进 FCM 算法 进 
行 图 像 分 割 。 首 先 ， 利 用 核 函 数 与 马 氏 距离 对 目标 函数 进行 
改进 ; 然后 ， 利 用 改进 后 的 算法 对 图 像 进行 分 割 。 通 过 对 比 
实验 证 明 ， 本 文 提 出 的 算法 对 于 人 工 合成 图 像 和 Berkeley 图 
像 库 图 像 都 能 有 效 解决 模糊 聚 类 算法 抗 品 性 问题 


1 FCMKM 算法 


1.1 核 函 数 与 马 氏 距离 
核 函 数 031 就 是 将 原 有 空间 的 像素 通过 新 的 特征 向 量 模 
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1 2 
p=-7> (3) 
过 2 
其 中 : c 为 聚 类 数目 ; n 为 样本 数 ;， 她 为 隶属 度 和 矩阵 U 的 元 
素 。 
b) Xie-Beni 系数 计算 公式 如 式 〈4) 所 示 。 
， > Fb 
HH, = 生 = 二 呈 性 - (4) 
名 -ul 
其 中 : I 表示 Xie-Beni 系数 的 紧 致 性 ,用 于 评价 类 内 的 内 聚 


式 映射 到 高 维特 征 空 间 ， 将 低 维 度 的 非 线 性 信息 转换 为 高 维 


度 线 性 问题 进行 处 理 。 对 所 有 的 x*，zEX， 函 数 K(x，z) 
=9(X)"9(z), 则 称 KK(x, z) 为 核 函数 , (x) 为 映射 函数 , p(x)-9(z) 
为 x、z 映 射 到 高 维特 征 空间 上 的 内 积 。 由 于 映射 函数 十 分 复 
杂 难 以 计算 ， 在 实际 中 通 带 者 站 使 用 核 函数 米 求解 内 积 计 
算 复杂 度 并 没有 增加 ， 映 射 函数 仅仅 作为 一 种 逻辑 映射 ， 表 
征 着 输入 空间 到 高 维特 征 空间 的 映射 关系 。 
常用 的 核 函 数 包括 线性 核 函数 、 多 项 式 核 函 数 、 高 斯 核 


度 ; Ze 表示 分 离 性 ， 用 于 评价 不 同类 间 的 耦合 度 。 因 为 好 的 
聚 类 划分 应 使 类 间 的 耦合 度 尽 可 能 低 ， 而 类 内 的 内 聚 度 尽 可 
能 高 。Xie-Beni 系数 越 小 ， 聚 类 的 内 聚 度 越 高 ， 耦 合 度 越 低 ， 
分 类 效果 越 明显 。 wj 为 隶属 度 秆 车 U 的 元 素 ; 、V、 了 从 为 聚 
类 中 心 矩阵 Y 的 元 素 ; 性 为 实验 对 象 的 元 素 ; m 为 模糊 系数 ; 
n 为 样本 数 ; c 为 聚 类 数 。 

c) 重 构 错误 率 用 来 刻画 分 割 后 的 重 构图 像 与 原 图 像 之 间 


函数 、Sigmoid 核 函 数 、 字 符 串 核 函数 。 本 文采 用 高 斯 核 函 
数 作为 实验 核 函数 ， 如 式 〈1) 所 示 。 


es 5] (1) 
其 中 ，w 表示 数据 的 标准 差 ， 其 对 数据 中 存在 的 噪声 具有 较 


好 的 抗 干扰 能 力 ， 但 却 具 有 很 强 的 局 部 性 。 为 改善 核 函数 局 
部 特性 ， 引 入 具有 全 局 性 的 马 氏 距离 ， 提 高 算法 的 全 局 性 ， 
避免 陷入 局 部 最 优 解 。 

马 氏 距离 0 (Mahalanobis distance ) 是 由 印度 统计 学 家 
马 哈 拉 诺 比 斯 (P. C. Mahalanobis) 提出 的 ， 表 示 数 据 的 协 方 
差距 离 。 它 是 一 种 能 有 效 地 计算 两 个 未 知 样本 集 相似 度 的 方 
法 。 与 欧 氏 距离 不 同 的 是 ， 它 通过 引入 协 方差 矩阵 ， 能 有 效 
地 描述 两 个 样本 点 之 间 的 全 局 性 关系 ， 包 含 更 多 的 邻 域 信息 
与 空间 信息 。 马 氏 距 离 计 算 公 式 如 式 〈2) 所 示 。 


M(x,y) = Ve -yD (x-y) (2) 


其 中 : 光 表 示 协 方差 矩阵 ， 为 x、y 之 间 的 协 方差 。 
1.2 聚 类 数目 与 初始 聚 类 中 心 选取 
聚 类 数目 和 初始 聚 类 中 心 的 选取 对 提高 模糊 聚 类 算法 的 
性 能 至 关 重 要 。 为 了 得 到 性 能 更 好 的 模糊 聚 类 算法 ， 本 文 对 
聚 类 数目 和 初始 聚 类 中 心 的 选取 方式 进行 改进 。 将 图 像 灰 度 
值 进行 映射 ， 得 到 灰 度 值 变化 折线 图 。 因 为 灰 度 值 变化 明显 
的 地 方 通常 对 应 目标 或 背景 所 在 的 区 域 ， 所 以 通过 计算 灰 度 
变化 折线 图 的 局 部 极 值 能 有 效 地 找到 最 佳 聚 类 中 心 。 将 局 部 
极 值 作为 聚 类 中 心 ， 局 部 极 值 的 个 数 作为 聚 类 数 。 具 体 的 计 
算 步 又 为 : 
a) 计 算 目 标 图 像 的 灰 度 直方 图 ， 
b) 计 算 灰 度 变化 折线 图 局 部 极 值 ， 
c) 将 局 部 极 值 作 为 聚 类 中 心 ， 
数 。 
1.3 评价 指标 
为 防止 主观 评价 因素 的 影响 ， 现 引入 Bezdek 划分 系数 
05、Xie_Beni 系数 09、 重 构 错 误 率 0 、 运 行 时 间 和 人 迭代 次 数 
五 个 指标 作为 评测 分 割 效果 的 指标 。 
a) Bezdek 划分 系数 是 一 个 判定 分 类 模糊 性 程度 的 标准 。 
分 类 越 分 明 时 ， 值 越 大 ; 分 类 越 模 糊 时 ， 值 越 小 。Bezdek 划 
分 系数 计算 公式 如 式 (3) 所 示 。 


得 到 灰 度 变化 折线 图 ， 
并 做 记录 ; 
局 部 极 值 的 个 数 作为 聚 类 


的 差别 ， 预 期 所 得 到 的 图 像 与 原 图 像 尽 可 能 地 相似 。 因 此 一 
个 好 的 图 像 分 割 算法 重 构 错误 率 应 尽 可 能 地 接近 0。 重 构 错 
误 率 计 算 公 式 如 式 (5) 所 示 。 


1 Dus 
R= 一 一 太 (5) 
i=l 2 
其 中 : zw 为 隶属 度 和 矩阵 U 的 元 素 ; c 为 聚 类 数目 ;xi 为 实验 
对 象 的 元 素 ; n 为 样本 数 。 


qd) 运行 时 间 是 指 算法 从 输入 图 像 到 图 像 分 割 完成 后 输出 
图 像 所 需要 的 全 部 时 间 。 通 过 计算 运行 时 间 ， 能 有 效 地 评价 
算法 的 时 间 复 杂 度 。 

e) 人 迭代 次 数 是 指 算法 在 更 新 隶属 度 和 聚 类 中 心 时 的 运行 
次 数 ， 每 更 新 一 次 迭代 次 数 加 1。 通 过 计算 迭代 次 数 ， 可 以 
有 效 地 评价 算法 的 收敛 性 。 
1.4 FCMKM 算法 
通过 对 模糊 聚 类 目标 函数 公式 〈 式 (6)) 进行 改进 ， 将 
核 函数 与 马 氏 距离 融入 式 (6) 中 ,通过 提高 模糊 聚 类 算法 邻 
域 信息 与 空间 信息 利用 率 ， 以 提高 算法 的 抗 噪 性 。 提 出 
FCMKM 算法 的 目标 函数 公式 ， 如 式 〈7) (8) 所 示 。 


Cou = 六》 | =—» 上 (6) 
本 1 拓 交 
Cc n 2 
G=> > ou MK sv,) (7) 


i=1 k=2 


MK (Gov) = KG, yD Kr) (8) 


其 中 : MK (xx，vi) 表示 第 个 元 素 到 第 i 个 聚 类 中 心 在 高 
维 空间 的 马 氏 距离 。 

根据 拉 格 朗 日 乘 数 法 , 求 得 隶属 度 迭 代 公式 如 式 (9) 所 
示 ， 聚 类 中 心 迭 代 公 式 如 式 (10〉 所 示 。 


ef MK, (9) 
下 


CE 
一 互 
V sy n 
>) 
j=l 


(10) 


.00070v1 


20181 


国 
国 


V 


chinaX 


王 燕 ， 等 : 基于 核 函 数 与 马 反 距离 的 FCM 图 像 分 割 算法 


FCMKM 算法 具体 流程 如 下 : 


a) 对 待 分 割 图 像 进行 初始 化 处 理 ， 通 过 本 文 改进 算法 ， 
天 数目 c 和 初始 聚 类 中 心 WO)， 设 定 迭 代 终 止 阔 值 = 
及 迭代 计数 器 b=0。 
居 式 (9 ) 计算 或 更 新 隶属 度 矩 阵 U(b), 据 公 式 (10) 
中 心 矩 阵 VB+7。 如 果 |lGCb)-GC+lll< = ， 则 和 迭代 


求 得 聚 


b ) 根 寺 


更 新 聚 类 


终止 ， 如 果 未 达到 终止 阔 值 但 是 达到 最 大 返 代 次 数 ， 则 友 代 
终止 ， 若 未 满足 迭代 终止 条 件 ， 则 和 迭代 计数 器 b 加 1， 并 重 


复 本 步骤 。 
) 迭代 终止 后 ， 根 据 式 (7) 计算 使 目标 函数 最 小 的 隶 
属 度 和 矩阵 U 和 聚 类 中 心 矩 阵 w, 即 为 最 佳 隶 属 度 和 矩阵 U 和 最 
中 心 矩 阵 V。 
于 任意 像素 点 习 ， 计 算 其 对 应 的 最 佳 隶 属 度 和 矩阵 
Ui;( 即 最 佳 


C 


佳 聚 类 


d 


Wij, 将 
度 值 。 
法 终 1 


隶属 度 和 矩阵 U 中 的 第 i 行 元 素 ) 中 元 素 的 最 大 值 
应 的 最 佳 聚 类 中 心 vj 的 灰 度 值 作为 像素 点 x 的 灰 


EE 复 此 过 程 ， 直 至 所 有 像素 点 灰 度 值 都 更 新 完毕 ， 算 


FE， 输出 分 割 后 的 图 像 。 


FCMKM 算法 流程 如 图 1 所 示 。 


本 
输入 待 分 割 图 像 


v 
初始 化 处 理 : 采用 改进 算法 ， 确 定 聚 类 
数目 c 和 初始 聚 类 中 心 V(0)， 设 置 兴 代 
终止 阔 值 及 迭代 计数 器 b=0 


v 
计算 或 更 新 隶属 度 算 阵 U(b)， 
并 更 新 聚 类 中 心 矩 阵 V(b+7) 


b=b+1 


达到 终止 阔 值 或 
达到 最 大 迭代 次 数 ? 


是 

v 

计算 求 出 最 佳 隶 属 度 和 矩阵 U 
和 最 佳 聚 类 中 心 矩 阵 V 


v 


根据 最 佳 隶 属 度 矩阵 UV 和 最 佳 从 类 中 
心 矩 阵 V， 更 新 图 像 像素 点 灰 度 值 


v 
输出 分 割 后 的 图 像 / 


v 
人 结束 ) 
图 1 FCMKM 算法 流程 
Fig. 1 FCMKM algorithm flow 


2 ”实验 结果 及 分 析 
为 了 验证 本 文 提 出 算法 的 有 效 性 ， 实 验 采 用 两 类 图 像 作 


为 实验 对 象 ， 分 别 是 人 工 合成 图 像 和 Berkeley 图 像 分 割 数 据 


库 图 像 。 同 时 对 传统 FCM、 基 于 核 函 数 的 FCM(KWEFLICM )、 
基于 马 氏 距离 的 FCM (FCM-M) 进行 分 割 性 能 比较 。 实 验 
环境 为 Intel R'Core TM'is 处 理 器 , CPU 主 频 为 2.6 GHz, 8 GB 
内 存 ，Windows 10 系统 的 PC 机 ， 采 用 MATLAB (R2016b) 


编译 语言 


实现 。 


2.1 参数 设置 
为 选取 模糊 加 权 指 数 m 的 最 佳 值 ， 现 对 m 进行 不 同 赋 
值 (m=1.5，1.6,.….，2.1)。 针 对 不 同 的 m 值 ,利用 FCMKM 


算法 对 区 


所 示 。 


2 (a) 进行 分 割 实验 ， 实 验 结果 如 图 2 (b) ~ (h) 
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(c)m=1.6 (d)m=1.7 


恬 从 侍者 


(e)m=1.8 


(fjm=1.9 


(gjm=2.0 (m=2.1 


妈 2 ”模糊 系数 实验 结果 


Fig.2 Experimental results of ambiguity coefficient 


针对 m 不 同 取 值 ， 对 实验 


疆 


一 口 


果 进 行 Bezdek 划分 系数 、 


Xie_Beni 系数 、 重 构 错误 率 三 个 客观 评价 指标 分 析 ， 结 果 如 


表 1 所 示 。 
表 1 模糊 系数 实验 结果 
Table 1 Experimental results of ambiguity coefficient 
Bezdek 划分 系数 Xie Beni 系数 重 构 错 误 率 

m=1.5 0.9988 0.8059 0 
m=1.6 0.9988 1.0496 0 
m=1.7 0.9989 0.5875 0 
m=1.8 0.9994 0.7805 0 
m=1.9 0.9996 0.8084 0 
m=2.0 0.9989 1.0573 0 
m=2.1 0.9991 1.0561 0 


根据 图 2 实验 效果 可 以 看 


H， 当 m=1.7 和 m=1.8 时 ， 算 
法 分 割 效果 较 差 , 因此 舍弃 m=1.7 和 m=1.8。 根 据 表 1 数据 


选取 Bezdek 划分 系数 较 大 、Xie Beni 系数 较 小 的 值 ， 综 合 
两 者 的 因素 ， 选 取 Bezdek 划分 系数 最 大 、Xie_Beni 系数 相 
对 较 小 的 m=1.9 作为 实验 用 模糊 加 权 指 数 。 

设置 迭代 终止 闪 值 := 104, 为 防止 迭代 多 次 后 仍 达 不 到 


和 迭代 终止 浆 值 的 情况 发 和 


2.2 人工 合成 图 像 


E， 设 置 最 大 迭代 次 数 为 100 次 。 


本 文采 用 图 3 中 四 幅 人 工 合成 灰 度 不 均匀 图 像 作为 实验 


图 像 。 


(a)8 字 医 


(b) 模 拟 噪 声 图 像 


(c) 球 形 图 案 (dd) 组 合 图 形 


妈 3 ”人 工 合成 图 像 


Fig.3 Synthetic image 
分 别 采用 传统 FCM 算法 、KWFLICM 算法 、FCM-M 算 


法 和 FCMKM 算法 四 种 算法 对 图 


结果 如 图 4 所 示 。 
通过 实验 结果 发 
FCM-M 算法 都 无 法 将 目标 和 


像 进行 分 割 实验 ,图 像 分 割 


了 岗 ， 传 统 FCM 算法 、KWFLICM 算法 、 


背景 有 效 地 进行 分 离 。 本 文 提 


少 了 53.99%， 是 


出 的 方法 能 有 效 地 将 目标 与 背景 


制 效果 提高 明显 ， 能 有 效 地 利 月 


分 离 ， 且 能 将 有 效 地 减少 噪 


声 的 干扰 。 为 防止 主观 评定 因素 的 干扰 , 利用 Bezdek 划分 系 
数 、Xie_Beni 系数 、 重 构 错误 率 三 个 指标 对 分 割 结果 进行 客 
观 分 析 ， 分 析 结 果 如 表 ~4 所 示 。 
通过 对 表 2~4 数据 


的 对 比 发 现 ,本 文 提出 的 FCMKM 算 
法 Bezdek 划分 系数 平均 提高 了 24.2%, Xie_Beni 系数 平均 减 
E 构 错误 率 平 均 减少 了 97.29% 。 本 文 算法 的 分 
邻 域 信息 与 空间 信息 ， 提 高 


算法 的 抗 噪 能 力 ， 对 噪声 图 像 分 割 效果 明显 。 


2.3 ” 抗 噪 性 能 实验 


为 检验 本 文 提 出 算法 的 抗 品 性 ， 对 图 3 (a) 进行 噪声 模 
拟 实验 。 对 图 3 (a) 分 别 添加 高 斯 噪声 、 椒 盐 噪 声 和 白 噪 声 ， 
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结果 如 图 5 所 示 。 
(a) 传 统 FCM  (b)KWFLICM (FCM-M (dJFCMKM 
1 最 a 《2 
的 【ms 
(e) 传 统 FCM KWFLICM (gDFCM-M (h)FCMKM 
© 
个 7 O90 
oee eee eee eee 
0 传统 FCM 0)KWFLICM (OFCM-M (DFCMKM 
和 
(m) 传 统 FCM nm)KWFLICM (oO)FCM-M (pFCMKM 
图 4 人 工 合成 图 像 分 割 效 果 对 比 
Fig.4 comparison of segmentation effect of artificial synthetic image 


表 2 Bezdek 划分 系数 实验 结果 


Table 2 ” Experimental results of Bezdek partition coefficient 


Eb IA 
能 ， 实 验 


算法 8 字 图 像 ”模拟 噪声 图 “球形 图 组合 图 形 
传统 FCM 0.8819 0.8926 0.8926 0.9247 
KWFLICM 0.9908 0.9987 0.9987 0.9997 

FCM-M 0.6583 0.9682 0.3438 -0.9977 
FCMKM 1 0.9941 1 0.9996 


表 3 Xie_Beni 系数 实验 结 


Table 3 Experimental results of XIe_ Beni coefficient 


算法 8 字 图 像 ” 模拟 噪声 图 “球形 图 组 合 图 形 
传统 FCM 1.2669 0.8993 0.8993 1.6522 
KWFLICM 0.7108 4.017 4.017 0.5604 

FCM-M 2.1564 0.4756 0.4958 -0.0322 
FCMKM 0.271 0.4284 0.1928 C0.0799 


表 4 重 构 错误 率 实验 结果 


Table 4 Experimental results of reconfiguration error rate 


算法 8 字 图 像 ”模拟 噪声 图 ”球形 图 组合 图 形 
传统 FCM 。 3.66E-28 1.49E-28 1.49E-28 1.91E-28 
KWFLICM 1.02E-29 2.84E-29 2.84E-29 1.53E-30 

FCM-M 1.10E-28 6.33E-29 6.81E-28 1.80E-28 
5.64E-30 0 0 


FCMKM 0 


图 5 模拟 噪声 实验 结果 


Fig.5 Experimental results of simulated noise 


: 基于 核 函数 与 马 色 距离 的 FCM 图 像 分 割 算 法 


将 添加 噪声 后 的 图 像 进行 分 割 ， 检 验算 法 的 抗 噪 性 


虹 声 后 的 图 像 ， 
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噪声、 椒盐 噪声 和 白 
(b)(g) (1) 为 普通 算法 对 三 种 噪声 图 像 分 


图 5 (a)(f)(k) 分 别 为 添加 高 #1 


果 。 


噪 性 
分 割 ， 且 背景 无 明显 噪声 ， 
验 表 明 ，FCMKM 算法 能 


割 结果 ; 
割 结果 ; 


结果 ; 


(c) (h) (m) 为 KWFLICM 算法 对 三 种 噪声 图 像 分 
Cd) (i) Cn) 为 FCM-M 算法 对 三 种 噪声 图 像 分 割 
(e) (j) (o) 为 FCMKM 算法 对 三 种 噪声 图 像 分 割 结 
通过 对 比 实验 结果 发 现 ，FCMKM 算法 能 有 效 地 提高 抗 
， 能 将 添加 高 斯 噪声 、 椒 盐 噪声 和 和 白 噪 声 后 的 目标 完整 
与 无 噪声 图 像 分 割 结果 相同 。 实 
了 效 地 提高 抗 噪 性 。 


BE 


2.4 ”Berkeley 数据 库 图 像 分 割 实验 
Berkeley 图 像 库 是 常用 的 图 像 分 割 数据 库 ， 利 用 该 库 中 
的 图 像 能 有 效 检 验 图 像 分 割 算法 的 性 能 。 为 验证 算法 分 割 效 
果 ， 选取 编 号 #159091、#3063、#15088 三 幅 图 像 作 为 实验 图 
像 ， 分 别 如 图 6 (a) ~ (c) 所 示 。 采 用 四 种 算法 分 别 对 三 幅 
图 像 进行 分 割 ， 实 验 结 果 如 图 6 (d) ~ (0o) 所 示 。 
(a) 能 原 
Lo 
Oh 传统 FCM WU EE M (g)JFCMKM 
由) 传统 FCM ed (OFCMKM 
(上 传统 FCM (m)KWFLICM nFCM-M (oO)FCMKM 


目标 


有 效 地 分 离 
景 进一步 地 分 离 ， 但 是 也 对 目标 造成 一 定 程度 的 模 煌 


图 6 Berkeley 图 像 实 验 结果 

Fig.6 Experimental results of Berkeley image 
通过 图 6 的 实验 结果 对 比 可 以 看 出 ， 传 统 FCM 算法 能 
标 分 离 出 来 ， 但 是 仍 保留 大 部 分 背景 ， 使 得 目标 与 背景 
; KWFLICM 和 FCM-M 算法 能 将 目标 与 背 
使 得 
完整 性 得 不 到 保留 ， 与 其 他 三 种 FCM 算法 相 比 ， 本 文 


提出 


的 FCMKM 算法 能 有 效 地 将 目标 与 背景 分 离 , 只 保留 
像 ， 使 得 分 割 结果 为 后 续 的 图 像 处 理 提供 方便 。 为 验证 


ej 
A 


的 有 效 性 ， 利 用 Bezdek 划分 系数 、Xie_Beni 系数 、 重 


数 平 
构 错 


FCMKM 算法 能 


为 错误 率 三 个 指标 对 图 像 分 割 结 
5 所 示 。 


进行 客观 分 析 ， 结 果 如 表 


通过 表 5 数据 的 对 比 发 现 , FCMKM 算法 Bezdek 划分 系 
均 提 高 了 48.3%，Xie_Beni 系数 平均 减少 了 53.5%， 重 
误 率 接近 0。 通 过 数据 分 析 可 以 看 出 ， 本 文 提 出 的 
有 效 地 进行 聚 类 划分 ， 能 提高 聚 类 内 部 的 内 


聚 度 ， 降 低 类 间 的 耦合 度 ， 降 低 重 构 错误 率 。 
2.5 FCMKM 算法 性 能 实验 


(b) 


代 次 
迭代 


为 验证 本 文 算法 的 分 割 效率 ， 对 图 3 (a) (b)、 图 6 (a) 
进行 运行 时 间 与 迭代 次 数 对 比 实验 ， 结 果 如 表 6 所 示 。 
从 表 6 可 以 看 出 ， 本 文 提出 的 算法 相 较 于 其 他 算法 在 迭 
数 与 运行 时 间 上 大 幅度 减 小 , 运行 时 间 平 均 减 少 77.3%， 
次 数 平均 减少 68.3%， 算 法 的 性 能 提高 明显 。 


录用 定稿 重 


燕 ， 等 : 基于 核 函 数 与 马 反 距离 的 FCM 图 像 分 割 算 法 
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表 5 Berkeley 图 像 分 割 评价 实验 结果 


Table 5 Experimental results of Berkeley image segmentation and evaluation 


能 图 像 


Bezdek 划分 系数 Xie_Beni 系数 重 构 错 误 率 Bezdek 划分 系数 Xie_Beni 系数 重 构 错误 率 Bezdek 划分 系数 Xie_ Beni 系数 重 构 错误 率 


传统 FCM 0.5719 0.6022 2.18E-28 0.5743 

KWFLICM 0.9932 1.0598 3.26E-29 0.9985 
FCM-M 0.8493 1.8855 1.75E-28 0.4073 
FCMKM 1 0.2816 0 1 


飞机 图 像 游船 图 像 
0.5359 3.20E-28 0.5288 0.0258 4.55E-28 
0.4912 3.99E-30 0.9786 0.4457 3.69E-29 
8.2309 1.84E-28 0.7508 0.8487 2.44E-28 
0.3074 0 0.9944 0.0489 0 


表 6 算法 性 能 实验 结果 实验 结果 


Table6 Experimental results of algorithm performance 


算法 图 3(a)8 字 图 。 图 3(b) 模拟 噪声 图 像 ”图 6(a)#159091 图 6(b)#3063 


迭代 次 数 运行 时 间 和 迭代 次 数 ”运行 时 间 和 迭代 次 数 运行 时 间 迭代 次 数 运行 时 间 
传统 FCM 65 6.971s 100 15.622s 19 15.047s 38 50.875s 


KWFLICM 4 1.135s 5 1338s 4 136s 100 200.072s 
FCMM 99 15.555s 100 15.045s 11 9.548s 53 77.079s 
FCMEM 3 1ll9s 5 0723s 8 5.783s 4 5.438s 
3 ”结束 语 


本 文 提出 了 一 种 结合 核 函数 与 马 氏 距离 的 模糊 聚 类 算 
法 ， 即 FCMKM 算法 。 通 过 对 图 像 进 行 核 函数 映射 ， 将 图 像 
像素 点 由 低 维 空间 非 线 性 映射 到 高 维 空间 ， 使 像素 点 变 得 线 
性 可 分 。 将 传统 的 欧 氏 距离 替换 为 马 氏 距离 ， 提 高 算法 对 邻 
域 信息 与 空间 信息 的 利用 率 。 通 过 对 比 实验 发 现 ， 本 文 提出 
的 算法 能 有 效 地 将 背景 与 目标 分 离 ， 提 高 了 模糊 聚 类 算法 的 
抗 噪 性 和 分 割 效 率 。 但 是 本 文 算法 也 存在 缺陷 ， 算 法 采用 单 
一 的 高 斯 核 函数 ,在 分 割 过 程 中 噪声 目标 部 分 细节 信息 丢失 。 
采用 自 适 应 的 核 函数 以 提高 算法 的 分 割 性 能 ， 是 下 一 步 的 研 
究 重 点 
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